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In recent years, the decay rj — > 3k has received considerable attention from experimental and 
theoretical side. It is of particular theoretical interest because it allows for the determination of 
(m u — mi). In addition, for many years now theory has had difficulties to understand the slope 
of the neutral channel Dalitz plot distribution, which in contrast is very well measured. I discuss 
here the relation of the decay to the masses of the light quarks and review a number of theoretical 
and experimental works that are concerned with these questions. 
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1. Introduction 

The masses of the light quarks, i.e., the up-, down-, and strange-quark, are not directly ac- 
cessible to experimental determination due to confinement, which prevents quarks from appearing 
as free particles. Since these masses are much smaller than the typical scale of QCD at around 1 
GeV, their contribution to hadronic quantities as, e.g., the nucleon mass is rather small. There is, 
however, a prominent exception to this rule, formulated in the famous Gell-Mann-Oakes-Renner 
relation [jl|]. It states that the masses of the lightest mesons are determined by the combined effects 
of spontaneous and explicit chiral symmetry breaking, that is by the chiral quark condensate and 
the light quark masses. At leading order in the quark mass expansion and including tj-tt mixing 
and first order electromagnetic corrections, one finds 

B {m u + m d )+A^ m + ... , 
B (m u + m s )+Af m + ..., (1.1) 



The parameter e = V3/4(m d — m u )/{m s — m) « 0.015, m = (m u + m d )/2, is the T7-7T mixing 
angle. According to Dashen's theorem [||], the electromagnetic corrections to the pion and kaon 
mass coincide at leading order with Ag/ n K ~ (35 MeV) 2 . Because Bo is not a priori known, one can 
only extract ratios of quark masses from Eq. (|TTJ). Assuming Dashen's theorem to be true, one 
finds the famous relations first derived by Weinberg using current algebra [||] : 

m d _ m 2 K0 -m 2 K+ +ml + m, _ m 2 K+ + m\ - m 2 n+ ^ 

Access to m u — m d is made difficult by the fact that this small quantity enters Eq. (JTTTJ) only quadrat- 
ically and by the possibility of Dashen violating contributions to the charged kaon mass. 

Lattice calculations are able to relate the quark masses to measurable meson masses, thus 
leading to reliable predictions for m s and m. I will not discuss these methods further and instead 
refer to the many contributions in these proceedings (e.g. as well as to the extensive report 

by FLAG 0. The determination of (m u — m d ) is however difficult also in Lattice calculations, due 
to the reasons discussed above: it enters quadratically or is hidden behind sizable electromagnetic 
corrections, which are only beginning to be studied on the Lattice [§, [J. 

2. 7] -^37T 

The main focus of this article is the decay process rj — > 3n. It can appear in two variants: 
either the decay goes into three neutral pions, Tj — > 3n°, or into a pair of charged pions together 
with a neutral one, tj — > 7t + 7t~7t°. I will denote the amplitude by £/ n (s,t,u) for the neutral and 
by srf c {s,t,u) for the charged channel. The Mandelstam variables are defined as s = {p K + + p K -), 
t = (p^ + p K -), and u = (p n o + p n +) in the charged channel, where they satisfy the relation s + 
t + u = m 2 1 + 2m 2 n+ + m 2 ^. Due to charge conjugation symmetry, the amplitude is symmetric under 



2e 



m^o = B (m u + m d ) + —B (m u - m d ) + 
|o =B (m d + m s ) + ... , 



2 n m u + m d + 4m s 2e 
in;, =B Q -j=B (m u - m d ) + . . . 
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t u. The adaptation of these definitions to the neutral channel is obvious. The neutral amplitude 
is even totally symmetric in s, t, and u. 

The particular importance of the decay v\ — y 37T for the determination of (m u — mj) is due 
to the fact that it is forbidden by isospin symmetry The reason is that three pions can not at the 
same time couple to vanishing angular momentum and zero isospin. The only operator in the QCD 
Lagrangian that can produce such a transition is 

Jz?/b = ^2 — -(uu — dd). (2.1) 

As a consequence of being generated by this AI = 1 operator, the decay amplitude must be pro- 
portional to [m u — m.d) and can be used to extract this quantity. The decay width can be seen as a 
measure for the size of isospin breaking in QCD. 

Due to the different electric charges of the up- and the down-quark, also the electromagnetic in- 
teraction is isospin violating and can contribute to the decay width of T] — > 37T. These contributions 
have been predicted to be small [|l(| 11], which has been confirmed by explicit one-loop calcula- 



tions within Chiral Perturbation Theory (ChPT) [|12J |13J. Let me stress, however, that Ref. Ql 3[ ] 
indicates that they might not be entirely negligible. Still, v\ — > 3tt provides a rather clean access to 
isospin breaking within QCD and hence to (m u — m^). 

In ChPT, quark masses are always multiplied by Bq. To avoid the appearance of this parameter, 
it is convenient to rewrite the prefactor to the amplitude in terms of quark mass ratios. Two different 
conventions are in use: 

« B (m u - m d ) = __L ^I«-^) + = m \) + , (2-2) 



Q 1 mi " ' R 



with 



2 ~2 

9 m—m m, — m 
Q 2 = — 2 j, R = — • (2-3) 



Note that in Eq. ( |2.2D the term containing Q is accurate up to G\M ), while the other one is 
corrected already at ^(^# 2 ). From these definitions, one finds that the two ratios are related by 



H2 V 



m 



-l 



R = 2Q 2 1 + • (2.4) 



The ratio m s /m can be determined on the lattice. The current lattice average from FLAG is 
m s /m = 27.4 ±0.4 [J7|]. Since the decay width is basically given by the phase space integral over 
the amplitude squared, 

f 2 1 1 

rV+3; t oe J \^^3n(s,t,u)\ — ^ , (2.5) 

an accurate theoretical description of the decay amplitude can be used to extract either Q or R by 
comparison with an experimental value for r^-^. Note that this also means that finding a value 
for Q or R is equivalent to finding the correct normalisation of the decay amplitude. 

Of course the aforementioned procedure can be reversed: given a reliable value for Q (or R) 
the decay width can be predicted. The theoretical challenge in doing this stays the same: one needs 
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an accurate description of the decay amplitude. The obvious choice to treat the process in question 
is of course ChPT. But this is not entirely successful, as one sees immediately by comparing the 
tree-level (or current algebra) [|4], [TJ), one-loop [16] and two-loop prediction [|7|] for the decay 



width 1 with the current PDG value, r^^+^o = 296 ± 16 eV [|T 

r^+^o =66eV + 94eV+138eV + ... = 298eV + ... . (2.6) 

Even though the result at two-loops happens to coincide almost perfectly with experiment, the 
chiral series does not exhibit good convergence behaviour. Since the decay width is based on 
the square of the amplitude, Eq. ( ^ ) is somewhat misleading: the numbers are enhanced by 
interference such that, e.g., the two-loop number also contains contributions that are of @{p % ) 
and @{p w ). The amplitude is converging more quickly. 

It has been shown independently of ChPT that the width is enhanced by large final state rescat- 
tering effects [|^]. This is mirrored in the chiral series and is a motivation to treat the process with 
dispersive methods, which allow to sum up the contributions from final state rescattering. 

Furthermore, the theoretical understanding of the energy distribution in the neutral channel 
is not complete. Due to symmetry reasons, the square of the amplitude can in this situation be 
parametrised by a single parameter a, which has been measured by many experiments with good 
mutual agreement. While experiments all find this parameter to be negative, ChPT at one- and two- 
loop order predicts a positive value. Also other calculations have failed to reproduce this parameter 
satisfactorily as can be seen in the compilation of results in Fig. ^. 

But difficulties appear also on the experimental side. There are hints of a tension between 
the many measurements of the neutral channel and the only available high-statistics measurement 
of the charged channel by KLOE. Such a statement can be made because the charged and neutral 
channel decay amplitudes are related by 

s/ n (s,t,u) = &f c (s,t,u) + &f c (t,u,s) + £/ c (u,s,t) , (2.7) 

if only first order isospin breaking, i.e. AI = 1, is considered. This relation allows for certain 
consistency checks among experiments, which will be discussed in more detail later. 

To conclude this section, I want to discuss an important property of the decay amplitude, the 



so-called Adler zero Q20|]. This soft-pion theorem states that in the St/ (2) chiral limit, the decay 



amplitude has two zeros at 

Pn+ -^0 ^ s = u = 0, t = m^ , and p n - -^0^s = t = 0,u = m^. (2.8) 

Moving away from the chiral limit (but keeping m n + = m^o = m n ), the positions of the Adler zeros 
are shifted by a contribution of the order of m\ to 

4 2 2 , ml . 2 ,._ ? , ml 



u = —m 



t = ml + — - , and s = t = -ml , u = ml + -^. (2.9) 



3 11 3 ' — 3 711 11 3 



'The tree-level and one-loop values have been taken from Q16|, where Q = 24.15 is used. This value follows from 
Dashen's theorem with the PDG meson masses from the time the article was published. No two-loop result for the decay 
width is quoted in JT^ ] since there, R has been calculated from the experimental value for T. But from their results 
together with Q = 24. 15 and m s /m = 21 A, one finds T = 298 eV. 
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As one expects from a SU (2) soft-pion theorem, the corrections to the position of the Adler zero 
are of order rrr % , since the symmetry forbids large @(m s ) contributions. At one-loop order, the real 
part has Adler zeros at s = u = 1.35m| and s = t = l.35m 2 , where also the imaginary part of the 
amplitude is small. 



3. Dalitz plot measurements 

The momentum distribution of a three-particle decay is typically displayed in form of a Dalitz 
plot, where it is plotted against two independent kinematic variables. A common choice are the 
so-called Dalitz plot variables, which in the charged channel are defined by 

X = 2^|>-"' ¥ = ^Q c « m "- m ^-^- 1 ' <31) 

with Q c =mjj — 2m n + — m^p . In the neutral channel it is for symmetry reasons convenient to use 
the variable Z = X 2 + Y 2 . The physical region of the decay process lies for both channels inside 
the circle with X 2 + Y 2 = 1 but does not entirely cover it. The point X = Y = is referred to as the 
centre of the Dalitz plot. An example of a three-dimensional Dalitz plot is shown in Fig. [T]. 

It is common to parametrise the Dalitz plot distribution as a polynomial in X and Y. For the 
charged channel and up to cubic order, the Dalitz plot parametrisation reads 

r c (X,Y) = \£/ c (s,t,u)\ oc 1 + aY + bY 2 + cX + dX 2 + eXY + fY 3 + gX 3 + hX 2 Y + IXY 2 , (3.2) 

where the coefficients a,b,... are called Dalitz plot parameters. Charge conjugation symmetry 
requires terms odd in X to vanish such that c = e = g = l = 0. 



Several experiments have measured at least some of the Dalitz plot parameters in Eq. (3.2) 



[22-24, 21], but only the measurement by KLOE has collected enough statistics to be reliable for 
more than one or two of them. Their result for the Dalitz plot distribution from about 1 .3 • 10 6 
v\ — > ti + ti~ti q events is shown in Fig. [I]. They find Dalitz plot parameters consistent with charge 
conjugation. Also the parameter h is consistent with zero. For the others, they find 

a = -1.090±^ , b = 0.124 ±0.012, d = 0.057^;^, / = 0.14±0.02. (3.3) 

This is the first time that / has been measured. 
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The Dalitz plot measurement of the charged channel by the WASA-at-COSY collaboration, 



which has been announced at this conference [|25||, has since been published in a PhD thesis [[2q] . 



The statistics is considerably smaller than for KLOE. The new data confirms the low value for b, 
but deviates from KLOE in d. Also the CLAS collaboration has collected large statistics on the 



charged channel, but the analysis has not yet been completed [27]. Furthermore, a new analysis of 



a much larger data set from KLOE is under way that also intends to improve on certain limitations 



of the previous analysis (see the contribution by L. Balkestahl in Ref. Q28|]). 

The neutral channel has been measured much more often in recent years [29-[37|] such that 
the picture is considerably clearer. Due to the symmetry in s, t, and u, the number of possible terms 
in the Dalitz plot parametrisation is much smaller in this case. Up to fourth order, it reads 

T n (X,Y) = \£/ n {s,t,u)\ 2 °c l+2aZ + 6/3y ^ 2 -y^ +2yZ 2 . (3.4) 

So far, experiments have not reached the accuracy needed to determine j8 and y, but a has been 
measured by many experiments in good agreement, leading to the average a = —0.0315 ±0.0015 
[18]. All the measurements entering this number are compiled in Fig. |J 

4. Theoretical work 

During the last few years, v\ — y 3ti has also received considerable attention from the theoretical 
side. I will in the following briefly discuss a few of the relevant works. An exhaustive discussion 
of the literature available on this subject is however not within the scope of this short article. 

• The strong contribution in the isospin limit has been calculated up to the two-loop level in 



ChPT Q17|]. The corrections to the one-loop result are found to be sizable (see also Eq. ( |2.6[ )). 
Unfortunately, a large number of low-energy constants enter, some of which are not well 
known. The size of the p 6 contribution can therefore not be estimated reliably at the present 
stage. 

The ChPT amplitude is used to extract a value for the quark mass ratio R = 41.3. Using 



Eq. (2.4) together with the FLAG average for m s /m, this leads to Q = 24.2. Furthermore, the 
result for the neutral channel slope parameter is positive, a = 0.013 ± 0.032, but due to the 
large uncertainty also encompasses negative values. 



Based on Eq. (2.7), an upper limit for a in terms of charged channel Dalitz plot parameters 
is derived: 

a<^(b + d-^cA . (4.1) 

It is this relation that can be used to check the consistency of charged and neutral channel ex- 
periments. The Dalitz plot parameters from KLOE and WASA-at-COSY both lead to an up- 
per limit that is negative but larger than the PDG average for a. While for KLOE the relation 
becomes almost an equality (a < —0.029), the upper limit is larger for the WASA-at-COSY 
result (a < —0.0036). The reason for this is mainly the larger value for d in the latter case. 
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The complete electromagnetic corrections up to order p A and e 2 m q have been calculated 



within ChPT in Ref. Q13|]. Compared to an earlier similar calculation [|12|], the terms of order 
e 2 (m u — rrid) have been added. Contrary to the assumption made in the older publication, 
these terms are of comparable size to other e 2 m q contributions and can not simply be ne- 
glected. However, the total electromagnetic contribution is still shown to be small; of the 
order of a few percent in the amplitude, less than one percent in the value of Q. 

The decay has been analysed within the framework of non-relativistic effective field the- 



ory (NREFT) up to two loops [38]. The method has before been successfully applied to 
K — )• 3 71 [ 39 - 4 1 ] and consists basically of an expansion in small pion three- momenta around 
the centre of the Dalitz plot. The calculation requires two main inputs. The low-energy con- 
stants appearing in the tree-level NREFT decay amplitude are determined by matching to the 
one-loop amplitude from ChPT at the centre of the Dalitz plot. Final state nn rescattering 
is included in the NREFT calculation. The additional low-energy constants that appear are 



fixed from two Roy equation analyses ([42, 43] and J44p). In this way, a representation of the 



shape of the Dalitz plot distribution in both channels is constructed. Since the overall nor- 
malisation of the amplitude can not be reliably determined within the NREFT framework, 
no value for the quark mass ratio Q is given. A particular advantage of the calculation is that 
isospin-breaking effects have been included. 

It is found that rescattering effects lead to sizable corrections to the Dalitz plot parameters. 
One- and two-loop contributions are in general of similar size, thus emphasising the partic- 
ular importance of rescattering effects in this process. Especially in the case of a, where the 
already positive tree-level value is further shifted in the positive direction by the one loop cor- 
rection, the two-loop contribution is large and negative. This leads to a = —0.025 ± 0.005, 
in marginal agreement with experiment. Regarding isospin breaking, it is found that the most 
sizable corrections to the Dalitz plot parameters are kinematic effects due to the fact that the 
charged and neutral pion masses are not the same. 



Using the NREFT method, it is possible to turn the upper limit in Eq. ( |4.1| ) into an equality. 
This is achieved by expressing the required correction in terms of the Dalitz plot parame- 
ter a. Since information on r\ — > 3n enters the NREFT amplitude through the Dalitz plot 
parameters, no input from one-loop ChPT is needed in this case. The Dalitz plot parame- 
ters from KLOE lead to a = —0.059 ± 0.007 in clear disagreement with the PDG average 



as well as KLOE's own value Q37[]. As the possible source of the problem, the parameter 
b is identified, since ^nreft = 0.308 > ^kloe = 0.124. However, the new measurement 
by WASA-at-COSY indicates that d rather than b might be responsible. Indeed, from their 
Dalitz plot parameters one finds a = —0.033 ±0.03, which encompasses the PDG average. 
Clearly, more data is needed in order to finally settle this matter. 



The process has been treated with a mostly analytical dispersive approach []45|], where two 
rescattering processes are taken into account. It has to be noted that chiral power counting 
is strictly followed such that the number of subtraction constants coincides with the number 
of low-energy constants at two-loop order which is six. Also, the dispersive result can be 
matched exactly to the two loop-result by an appropriate choice of subtraction constants. 
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Figure 2: The one-loop amplitude from ChPT together with the dispersive amplitude from Ref. [ p3| | along 
the line s = u. Dashed lines represent real parts, dotted lines imaginary parts. It is clearly visible that the 
dispersive amplitude has no Adler zero. I thank K. Kampf for providing the program that produced this data. 

The main result of this work is a dispersive representation where the subtraction constants 
are fitted to the charged channel data from the KLOE collaboration. Since the overall nor- 
malisation depends on Q, it can not be obtained from the data. Instead, the imaginary part of 
the dispersive amplitude along the line t = u from zero up to threshold is required to be close 
to the two-loop result. The motivation for this procedure is that the 0(p 6 ) corrections happen 
to be rather small along this line. The resulting amplitude reproduces the experimental Dalitz 
plot distribution in the physical region and leads to Q = 23.1 ±0.7 and a = —0.0044 ±0.004. 

There is however a severe problem with the dispersive amplitude. As can be seen in Fig. ||, 
it has no Adler zero, which means that it is not consistent with SU (2) chiral symmetry. It is 
hard to justify any use of ChPT in a calculation that so blatantly violates the symmetry at its 
foundation. 



5. Our dispersive analysis 

I want to conclude the discussion of the theoretical literature on v\ — y 2>7i by a somewhat 
more detailed description of an ongoing numerical dispersive calculation in collaboration with 
G. Colangelo, H. Leutwyler, and E. Passemar. Intermediate status reports of this work have been 



presented at conferences before [ ]28| , f!6j p7| | and have been published in form of a PhD thesis [[T8|] . 
A more detailed discussion of the formalism can be found in these references and I will focus here 
on newer developments that are not described there. 



The method we apply has been described in detail in Refs. [ }49| , |50p. Because the calculation 
is numerical, it is possible to include an arbitrary number of rescattering processes. Also, chiral 
power counting is not followed: we always use the best available input, e.g., we do not expand the 
7171 scattering input to the appropriate order in each iteration. The method involves two main steps 
that can be treated entirely independently: A dispersive representation must be derived and solved 
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numerically. The result is an amplitude that is a function of a number of unknown subtraction 
constants. The second step is then to determine these constants in a good way. 

The dispersion relations are based on the decomposition of the normalised decay amplitude 

2 

J?(s,t,u) =M Q (s) + (s-u)M l (t) + (s-t)M 1 (u)+M 2 (t)+M 2 (u)- -M 2 {s), (5.1) 

where the subscript stands for the isospin of the scattered pion pair. This decomposition has the 
advantage that one deals with functions of a single variable only. Using analyticity and unitarity, 
one arrives at a set of dispersion relations for the three functions Mj(s) that are all of the form 

, . f , , s ni r ds' sin 5/(y)M/(j') ] 



where Q/(s) are the so-called Omnes functions Q52| ] given by 

= exp < 



' /./-' . * . > . (5.3) 



7T J 5' (5' — 5 — IE) 

Ami 



The functions are angular averages over all the Mj(s), such that the dispersion relations are 

recursive and coupled with each other. Two kinds of inputs are needed. On the one hand, one needs 



to know the nn scattering phase shifts 8j (s), which we take from Ref. Q43|]. On the other hand, the 
polynomials Pi(s) contain the subtraction constants, which must be determined with information 
from outside the dispersive machinery. It is on this issue that progress has taken place recently. 



The dispersive representations in Eq. ( p.2| ) can be Taylor expanded: 



Mi{s) = ai + bis + as 1 + djs 3 + ... . (5.4) 

The Mj(s) are uniquely determined once all the subtraction constants are fixed. They have unique 
Taylor expansions, such that one can define a unique relation among a set of subtraction constants 
and an appropriate equally sized set of Taylor coefficients. A solution of the dispersion relations can 
therefore be specified by giving values for the subtraction constants or for the Taylor coefficients. 
The advantage of working with the latter is that they pick up imaginary parts only at ^(p 6 ), such 
that they can safely be approximated as real. This then automatically leads to subtraction constants, 
where the imaginary part is non-zero, but suppressed compared to the real part. 

The splitting of the amplitude into the Mj(s) is not unique because of the relation s + 1 + u = 
+ 2m 2 K+ + m^ : one can add polynomials to theM/(j) in such a way that M(s,t,u) is not changed. 
This gauge freedom allows to choose some Taylor coefficients arbitrarily and the number of free 
parameters is therefore smaller than the number of Taylor coefficients (or subtraction constants) 
that are used. 

We have tried to work with different numbers of subtraction constants, but comparison with 
data has shown that satisfactory agreement is only achieved, if eleven constants are used. The 
corresponding Taylor expansions read 

M (s) =a + b s + c s 2 + d s 3 + . . . , 

M l (s)=a l +b l s+ s 2 + ..., (5.5) 
M 2 {s) =a 2 + b 2 s + c 2 s 2 + d 2 s 3 + ... . 
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The reason for ending the series one order lower for M\ (s) is that this function is multiplied by 
another power of momenta in Eq. (5T), such that all three functions lead to contributions of &(s ) 



in the total amplitude. 

We determined the subtraction constants with two different methods. The first one relies en- 
tirely on one-loop ChPT as input and makes no use of data at all: 

• Since one-loop ChPT can not be trusted for terms beyond &(s 2 ), the coefficients do, c\, and 
do are set to zero. The dispersion relations thus contain eight subtraction constants. 

• The parameters ao, bo, a\, and 02 are set to their tree-level value using the gauge freedom. 
This means that there are actually only four free parameters in the dispersion relations. 

• Finally, the remaining parameters CQ,b\,b%, and ci are set to their one-loop value. 

In this way, the dispersive solution is entirely fixed. It will later be referred to as "dispersive, one 
loop". I stress again that for this solution, only ChPT at low energy is used in order to fix four 
subtraction constants. This low-energy information is then extrapolated to the physical region by 
means of the dispersion relations. 

The second method uses the full set of eleven subtraction constants: 

• The parameters ao, bo, a\, ci2, Co, b\, bj, and C2 are determined exactly as before. 

• The presence of more parameters also enlarges the gauge freedom since polynomials of 
higher order can now be added to the Mj(s). The parameter di is chosen such that ^2, which 
is the coefficient to s 3 in Pi(s), is zero. The number of free parameters is thus six in this case. 

• Finally, the remaining two parameters and are determined by fitting the square of the 
amplitude in the physical region to the charged channel data from the KLOE collaboration. 

This leads to another solution of the dispersion relations which will be referred to as "dispersive, 
fit to KLOE". We have also fitted other available data sets, but for simplicity, only one of these fits 
is presented here as an example. 

The Dalitz plot distributions from both solutions are depicted in Fig. |3| together with the cor- 
responding one-loop results. The left panel shows the Y distribution along the line with X = 0, 
which is equivalent to the s distribution for t = u. Note that small s correspond to large Y due to 
the minus sign in the definition of Y in Eq. (JO|). Clearly, ChPT is successful at low energy, but 
fails towards the upper end of the physical region. It is noteworthy that the extrapolation of the 
low-energy ChPT amplitude through the dispersion relation already leads to a considerably better 
agreement with experiment, which is then further improved by the fit. It seems not, at first sight, 
remarkable that the fit does agree with the data it is fitted to. However, the fit must at the same 
time be consistent with the chiral constraints at low energy and it is important to show that this can 
actually be achieved. 

Through Eq. ( |2.7[ ), the neutral channel amplitude can be calculated from the charged channel. 
The right panel of Fig. |3| shows the Z distributions that one finds in this way from the three charged 
channel amplitudes. Accordingly, the blue curve involves no experimental information on the 
neutral channel. But the influence of the charged channel data is exactly what is needed to bring 
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Figure 3: Results from the numerical dispersive analysis. Left panel: The square of the amplitude in the 
charged channel along the line X = in comparison with the KLOE curve. The dashed line marks the end 
of the physical region. Right panel: The square of the amplitude in the neutral channel integrated along 
circles of constant Z in comparison with data points from MAMI-C. Note that none of the theoretical curves 
involves experimental information on the neutral channel. 



the shape of the Z distribution into agreement with the neutral channel data. To visualise this, the 
figure contains data points from the MAMI-C collaboration [|36|] as an example. 

But not all problems are solved yet. The dispersive amplitude is evaluated in the isospin limit, 
while the data are collected in the real world with two different pion masses and electromagnetic 
effects. These effects are expected to be small but if they are entirely neglected, we predict a neutral 
channel Dalitz plot distribution that is not in agreement with the data. We expect the largest isospin 
correction to be kinematic effects due to the pion masses. Indeed, taking these into account by 
shifting the corresponding singularities to their physical position through a slight deformation of 
the phase space, we find that a fit to charged channel data leads to good agreement with neutral 
channel data as well. 

While isospin corrections to the shape of the Dalitz plot distribution are successfully treated 
in this way, the decay rate suffers from the procedure. This can be seen in the fact that we find 
a branching ratio r = r^^^o /r^^^+^-^o that is not in agreement with experiment. But from the 
estimates for electromagnetic effects in Ref. [13 J, we expect that the one-loop result from ChPT 
including isospin breaking can be used to remedy the situation. 



6. Comparison of results for a and Q 

To conclude, I have compiled various experimental and theoretical results for a and Q in 
Figs. f| and |5[ In particular, all the results from works that I have mentioned above are listed. 

Figure ^| clearly shows the failure of ChPT at one- and two-loop order to reproduce even the 
sign of a. On the other hand, all four dispersive results do reproduce the sign, and our fit to KLOE 
is even in agreement with the PDG average. Also the NREFT calculation leads to a value that is 
compatible with experiment. The current PDG average, which is marked by the grey band, includes 
all the experimental values that are given in the figure. 

The values for Q that are shown in Fig. || cover the range from about 21 up to Qd = 24.3, 
which follows from Dashen's theorem. Our preliminary value from the fit to KLOE lies around 22. 
I have conservatively assumed an error of one unit, but expect the final error to be smaller. 
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Figure 4: Comparison of various theoretical and experimental results for the slope parameter a. 
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Figure 5: Comparison of various theoretical and experimental results for the quark mass ratio Q. 



7. Conclusion & Outlook 

The process T] — > 3n is a unique source of information about isospin breaking in QCD and thus 
for estimating m u — m^. A wealth of theoretical and experimental work has been dedicated to this 
decay and the question of quark mass extraction. With dispersive techniques, the sizable final state 
rescattering effects can be treated properly and a significant improvement of the one-loop result 
from ChPT has been achieved in this way. In particular, the neutral channel slope parameter can be 
understood based on charged channel data. Consequently, the dispersive analysis does not produce 
a clear sign of a tension among experiments, even though that possibility can not yet be entirely 
excluded. Before a final statement can be made, a more careful treatment of isospin breaking 
effects must be implemented. 
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